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The  inverses  of  r-banded  matrices,  for  r  =  1,2,3  have  been  thoroughly  investigated  as  one 
can  see  from  the  references  we  provide.  Let  B,  n  (1  <  r  <  n)  be  an  n  x  n  matrix  of  entries 
{ aj } ,  —r<i<r,l<j<r,  with  the  remaining  un-indexed  entries  all  zeros.  In  this  paper, 
generalizing  a  method  of  Mallik  (1999)  [5],  we  give  the  LU  factorization  and  the  inverse 
of  the  matrix  Br„  (if  it  exists).  Our  results  are  valid  for  an  arbitrary  square  matrix  (taking 
r  —  n),  and  so,  we  will  give  a  new  approach  for  computing  the  inverse  of  an  invertible 
square  matrix.  Our  method  is  based  on  Hessenberg  submatrices  associated  to  Br  n. 

©  2012  Elsevier  B.V.  All  rights  reserved. 


1.  Background 


There  are  various  methods  for  finding  inverses  of  matrices  (if  these  exist),  and  we  recall  the  Gauss-Jordan  method,  the 
triangular  decomposition  such  as  LUD  or  Cholesky  factorization,  to  mention  only  a  few.  A  very  popular  approach  is  based 

on  block  partitioning.  Let  A  =  ^  whose  inverse  (called  Schur’s  complement)  is 


A-1 


+  ^n^\2B  — A111A12B  1 

-fr'A^A-1  B-1 


where  B  =  A22  —  A2JAn  A12  (under  the  “right”  conditions,  namely,  that  An ,  B  are  invertible).  Further,  an  iterative  approach 
for  finding  the  inverse  exists,  but  its  convergence  is  a  function  of  the  matrix’s  condition  number,  cond(A)  =  ||Aj|  || A- 1 1|  (in 
some  norm,  say,  £oo,  where  HAIloo  =  max  ^  |a,j|;  or  £i,  with  ||A||i  =  || AT  || oo ).  We  will  not  go  into  details  regarding  these 
parameters  of  a  matrix  as  the  considerable  literature  has  been  dedicated  to  these  concepts. 

Certainly,  the  inverse  of  a  diagonal  matrix  (if  it  exists)  is  found  easily:  if  A  =  diag(a\, . . . ,  an),  then  A-1  = 

u-[  C\  0  •  •  •  0 


diag(a1 


,  an ').  If  U  is  a  bidiagonal,  say  U  = 


satisfy  the  recurrences 


a  2  c2 


,  Uj  r/z  0,  then  U 


(vij)ij,  where  the  entries 


0, 

1 


Ui’ 

Qfi+lj 

Ui 


i  >  j; 
i  =J 
i  <  j. 
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These  recurrences  can  certainly  be  solved,  and  we  obtain 


i  >j 
>  <j- 


Going  further  to  tridiagonal  matrices,  things  start  to  change.  We  want  to  mention  here  the  work  of  Schlegel  [  1  ],  who  showed 

~bi  1  0  •••  O' 


that  if  B  ■■ 


1 

0 


0 

bn 


,  then  the  inverse  B  1  =  is  given  by  eg  =  r  1  r,_  i  rn_i,  i  <  j  and  c,-,-,  if  j  <  i,  where 


r0  =  1,  ri  =  —bi,  n,  =  —(bijk-i  +  r&_ 2 ),  k  =  2, _ n  —  1  and  r  =  b„r„_i  +  rn_2  =  (— l)n+1  det(B).  We  would  like  to 

make  an  observation  at  this  point:  since  the  inverse  of  a  tridiagonal  matrix  is  a  full  matrix,  the  Schur’s  complement  method 
is  not  very  efficient. 

Moreover,  Vimuri  [2]  obtained  the  inverse  of  another  particular  tridiagonal  matrix  in  terms  of  the  Gegenbauer  polynomial 
C“  (x)  (for  a  —  1)  whose  generating  function  is  (1_2J+t2)1,  =  Y^Lo  (x)tn;  Prabhakar  et  al.  [3]  showed  some  connection 

between  the  aforementioned  inverse  and  the  generalized  Hermite  polynomials,  namely,  g™(x,  A)  =  k!(n-mk)!^kxn~mk' 

(Classical  Hermite  H„(x)  are  obtained  for  m  =  2,  X  =  —  1,  x  ->  2x.) 

Many  special  cases  of  the  banded  matrices  such  as  Toeplitz  matrices,  symmetric  Toeplitz  matrices,  especially  tridiagonal 
matrices,  etc.,  have  been  studied  by  several  authors,  as  we  previously  mentioned  (the  reader  can  find  many  more  references 
published  in  the  present  journal,  for  instance  [4],  among  others).  These  matrices  arise  in  many  areas  of  mathematics  and  its 
applications.  Tridiagonal,  or  more  general,  banded  matrices  are  used  in  telecommunication  system  analysis,  finite  difference 
methods  for  solving  PDEs,  linear  recurrence  systems  with  non-constant  coefficients,  etc.,  so,  it  is  natural  to  ask  the  question 
of  whether  one  can  obtain  some  results  about  the  inverse  of  4-diagonal,  or  perhaps,  even  general  banded  matrices.  We  will 
do  so  in  this  paper. 

Throughout  this  paper,  we  consider  a  general  r-banded  matrix  B,„  of  order  n  defined  by 


B,,n  =  ( by ) 


0  Cr+1 


0 

a. 


an-1 


n~  2 
un- 2 


n— r+1 


n— 3 
2 


(i) 


where  aj’s  stand  for  arbitrary  real  numbers,  —r<i<r,\<j<r<n  (note  that  the  i  in  the  notation  a'n  does  not  denote 
the  ith  power  of  an ).  When  r  =  n,  the  matrix  Bn  n  is  reduced  to  an  arbitrary  square  matrix. 

In  this  paper,  generalizing  a  method  of  [5  ],  we  give  the  LU  factorization  and  (in  our  main  result)  the  inverse  of  the  matrix 
Br  n  (if  it  exists).  Our  results  are  valid  for  an  arbitrary  square  matrix  (by  taking  r  =  n).  Therefore  we  give  a  new  approach  for 
computing  the  inverse  of  an  invertible  square  matrix,  thus,  generalizing  various  results  (see  [6-8,5,9,10,1  ],  and  the  references 
therein).  Our  method  is  based  on  Hessenberg  submatrices  associated  to  B,  „. 

Section  2  deals  with  the  LU  factorization  of  an  r-banded  matrix.  To  find  the  inverse  of  such  a  matrix  (obtained  in  Section  4), 
we  need  to  find  inverses  of  the  obtained  triangular  matrices  from  the  LU  factorization  and  this  is  done  in  Section  3.  We 
conclude  the  paper  with  a  few  examples  in  Section  5. 


2.  LU  factorization  of  an  r-banded  matrix 

This  section  is  mainly  devoted  to  the  LU  factorization  of  the  matrix  Br  n.  First,  we  construct  two  recurrences.  For  1  <  i  <  r 
and  s  >  r  >  1,  define 


-E> 


(2) 
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and  for  1  <  i  <  r  —  1 


a, 


— O+l) 


ml  = 


t=  i 


i+t  ht+ 1 
s—t 


k} 


(3) 


with  initial  conditions 
r-0+l) 

'1  _ 


m  = 


-0+D 

l 


k] 


1  <  i  <  r  -  1 


and 


1  <  i  <  r, 


where  the  terms  a^1  for  1  <  i  <  n  are  the  entries  of  Brn.  For  these  sequences  to  be  well-defined,  we  assume  that  none  of 
the  denominators  k's  are  zero  (which  is  equivalent  to  the  below-defined  U,  and  consequently,  Br  „  being  invertible). 

Now  define  unitary  lower  tridiagonal  matrix  L  and  upper  triangular  matrix  U, 


L  —  (lij)  — 


and 


1 

m? 


m 


m 


m. 


m\ 


mA 


n— r+l 


0- 


"h-2 

„2 


(4) 


U  =  ( uu )  = 


0 


k\ 


k22 


kl 


]A 

K  4 


l<\  0 


'<  4 


kr2 


kr„ 


n—2 

l  ,1 


n— r+l 


vn— 2 

k2„-t 

iA 


(5) 


Our  first  result  gives  the  LU  factorization  of  Br  , 


Theorem  1.  For  n  >  1,  the  LU  factorization  of  matrix  Br  jI  is  given  by 
B,,n  =  LU 

where  L  and  U  are  defined  as  in  (4)  and  (5),  respectively. 

Proof.  First,  we  consider  the  case  1  <i=j<r.  From  matrix  multiplication  and  the  definitions  of  L  and  U,  we  have 

n  i 

h;;  =  ^  ^  ^  ^  li.li., 

s=  1  s=  1 

=  klul.i  +  k2u2,i  +  •  •  •  +  liiUjj 

=  m[-'k\  +  m'-2/^1  +---  +  ml,kf_1  +  fc?. 

From  (2),  by  taking  r  =  i,  we  obtain 
bu  =  aj, 
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which  gives  the  conclusion.  Now  consider  the  case  of  i  >  r.  Thus 

n  r 

^  ^  hstisi  —  ^  s+l^i— s+\,i 

s=  1  s=  1 

=  IffUfi  +  — l Wf — iTf  +  *  •  •  +  ^i,i— r+l^i— r+l,i* 

From  the  definitions  of  L  and  U,  we  write 

b„  =  k,1  +  +  ml2kl2  +  ■  ■  •  +  m^k^ 

which,  by  taking  i  =  1  in  (2)  implies 
bu  =  a'n, 

which  shows  the  claim  for  i  =  j. 

Next,  we  look  at  the  super-diagonal  entries  of  matrix  Br  n.  Consider  the  case  j  =  i  +  q  where  1  <  q  <  r- 
definition  of  Bn,  for  1  <q  <r-  1, 


bi,i+q  —  ^  ^  h,sUs,i+q- 

s=  1 

We  consider  two  cases.  First,  we  assume  1  <  i  <  r  —  q,  and  so, 

i 

ktut,i+q  —  klUl,i+q  +  h,2u2,i+q  +  *  '  *  +  hiui,i+q 

t=  1 

=  +  ni2+9/4+q_1  H - F  +  /<•“ 

which,  by  taking  i  —>  q  +  1  and  n  — »•  i  in  (2),  gives  us 


J2ktUt.i+q  =  a-+q, 

t=  1 

which  completes  the  proof  for  the  first  case. 

Now  we  consider  the  case  r  —  q  <  i.  Thus  we  write 


bi,i+q  —  ^  ^  k,sUs,i+q  —  ^  ^  h,i—t 


t+l,i+q. 


From  the  definitions  of  matrices  U  and  L,  we  can  write 


r-q 

El‘ 

t=i 

,i—  t+l^i—  t+l,i+q  —  liiUi,i+q  H"  *i,i—  1  Wj—  l,i+q  "T  ’  *  *  H"  /j,; 

i-r-q+lui-r+q+l,i+q 

=  kf+q  +  m\_.,kffq  +  mf_2kf+q  +  ■ 

.  .  _i_  m'-9-1  // 

'  r+q+l'S— r+q+1 

=  kl+q  +  r)Tl'mitkqfw. 

t=i 


Using  (2),  we  obtain 

.  1 +q 

bi,i+q  —  Qj  , 

which  completes  the  proof  for  the  upper  diagonal  entries  of  matrix  Br 

Finally,  we  look  at  the  upper  diagonal  entries  of  the  matrix  Br  n,  and  we  need  to  show  that  bi+q  i  =  aj~(9+1) 
consider  the  case  1  <  i  <  r  —  q.  From  the  definitions  of  matrices  U  and  L, 


bi+q,i  —  ^  '  h+q,tUt,i  —  ^  '  h+q.t^tj 


=  m\+q  V,  +  mi2q  *k2 1  H - F  +  mqk, 


t=  i 
i+q— 2 ,  —  1 


1.  Thus,  by  the 


.  Flere,  we  first 
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which,  by  taking  n  ->  i  and  i  ->  q  in  (2),  implies 


h  ■  —  n~^q+^ 

Di+q,i  —  ai 

For  the  final  case,  that  is,  i  >  r  —  q,  using  the  definitions  of  U i  and  Llf  we  write 

n  r—q 

bi+q,i  =  ^  '  k+q,tUt,i  =  ^  '  h+q,i— t+jUj— t+l,i 
t=l  t=l 

=  k+qjMiJ  H"  li+qj—  \  l/j—  i;j  “h  *  *  *  “h  h+q,i— r+q+1  r+q+1,/ 

-  mqkx  +  mq+V  + _ h  m?+(r-<?-i),r-g 

—  fiijKj  -f  rn1_1  T-  T-  "*i_r+9+1  N— r+q+1 


t=i 


which,  by  taking  n  — >  i  and  i  — >■  q  in  (2),  gives 

and  the  theorem  is  proved.  □ 

The  result  of  Theorem  1  will  be  valid  for  the  LU  factorization  of  any  arbitrary  square  matrix  by  taking  r  =  n  in  the  matrix 
Now  we  give  a  closed  formula  for  det  ( Br  n )  using  the  LU  factorization  of  Br  n. 


Corollary  2.  For  n  >  0, 

n 

detBr,„  =  Y\ 

!=1 

where  l<l ’s  are  given  by  (2). 


3.  The  inverse  of  triangular  matrices 


In  this  section  we  give  an  explicit  formula  for  the  inverse  of  a  general  triangular  matrix.  For  this  purpose,  we  construct 
certain  submatrices  of  a  triangular  matrix,  which  are  Hessenberg  matrices,  and  then  we  consider  the  determinants  of 
these  submatrices  to  determine  the  entries  of  the  inverse  of  the  considered  triangular  matrix.  Since  upper  and  lower 
triangular  matrices  have,  essentially,  the  same  properties,  first  we  consider  the  upper  triangular  matrix  case.  We  denote 
the  corresponding  Hessenberg  matrices  for  an  upper  and  a  lower  triangular  matrix  by  Hu  (r,  s )  and  He  (r,  s),  respectively. 
Let  H  =  ( hij )  be  an  arbitrary  (n  x  n)  upper  triangular  matrix.  Now  we  construct  square  Hessenberg  submatrices  of  H 

of  order  |s  —  r|  in  the  following  way:  for  s  >  r  >  0,  let  Hu  (r,  s)  =  (  h jj)  denote  an  upper  Hessenberg  submatrix  of  H  by 

deleting  its  first  r  and  last  (n  —  s)  columns,  and,  first  (r  —  1)  and  last  (n  —  s  +  1)  rows.  Clearly  the  (s  —  r)  x  (s  —  r)  upper 
Hessenberg  matrix  Hu  (r,  s)  takes  the  form: 


hr,r+ 1  hr, r+2 

hr+  l,r+l  hr+  i,r+2 


Hu  (r,  s)  = 


0 


hr,s—  1  hrs 
hr+l,s—1  hr+ 1,J 


0 


hs-2,s-2 

0 


hs—2,s—  1 
h$-  i.s—  i 


hs-2,s 

hs-u- 


(6) 


Similarly  let  H  =  (fi,j)  be  an  arbitrary  (n  x  n)  lower  triangular  matrix.  We  construct  square  Hessenberg  submatrices  of 

H  of  order  (r  —  s)  in  the  following  way:  for  r  >  s  >  0,  let  Ht  (r,  s)  =  denote  a  lower  Hessenberg  submatrix  of  H  by 

deleting  its  first  r  and  last  (n  —  s)  rows,  and,  first  (r  —  1)  and  last  (n  —  s  +  1)  columns.  Clearly  the  (r  —  s)  x  (r  —  s)  lower 
Hessenberg  matrix  H?  (r,  s )  takes  the  form: 


'r+l.r 

ftr+l,r+l 

0 

0 

'r+2  ,r 

hr+2,r+\ 

hs-2,s-2 

!s-l,r 

hs—l,r+ 1 

hs-l,s-2 

fts—  l,s—  1 

hs,r 

hs,r+ 1 

hs,s—2 

^s,s— 1 
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Here  we  note  that  (Hu  (r,  s))T  =  (He  ( s ,  r)).  Our  construction  is  valid  for  any  upper  or  lower  triangular  matrix,  but  we  will 
apply  it  only  to  L,  U,  given  by  (4)  and  (5),  rendering  Le(i,j),  Uu(i,j). 

For  example,  let  A  be  an  upper  triangular  matrix  of  order  6  as  follows: 


_ai 

bi 

Cl 

d, 

Cl 

/n 

0 

a2 

h2 

c2 

d2 

?2 

A  = 

0 

0 

0 

0 

a3 

0 

b3 

a4 

c3 

fa4 

d3 

c4 

0 

0 

0 

0 

as 

fas 

_0 

0 

0 

0 

0 

a6- 

ThusAu  (2,  5)  and  Au  (3,  6)  take  the  forms : 


"fa2 

c2 

d2l 

'fas 

c3 

d3 

(2,  5)  = 

a3 

_0 

fas 

a4 

1 

ro  ^ 
O 

and  /4U  (3,  6)  = 

a4 

_0 

fa4 

a5 

c4 

fas. 

Here  note  that  all  matrices  of  the  form  Hu  ( i,j )  (Ht  ( i,j ))  obtained  from  an  upper  (lower)  triangular  matrix  U  are  upper 
(lower)  Hessenberg  matrices. 

Now  we  start  with  the  following  two  lemmas.  Throughout  this  paper,  we  assume  the  boundary  conditions  Hu  (r,  r)  = 
He  (r,  r)  =  1  and  x,  =  1  for  i  >  j. 


Lemma  3.  Let  the  (j  —  i)  x  (j  —  i)  upper  Hessenberg  matrix  Hu  ( i ,  j )  be  defined  as  in  (6).  Then,  for  j  >  i  +  1, 

{J-i 
(-l)2j-t  det  (Hu  ( i,j  -  t))  ]~[  akk 

k=j-\-t 

Proof.  If  we  compute  the  determinant  of  the  upper  Hessenberg  matrix  H„  ( i,j )  by  the  Laplace  expansion  of  a  determinant 
with  respect  to  the  last  column,  then  the  proof  follows.  □ 

Theorem  4.  Let  U  =  (a,j)  be  an  ( n  x  n)  arbitrary  upper  triangular  matrix  and  W  =  (wy)  =  U~\  its  inverse.  Then 
(a,-,)-1  ^  if  i  =  j, 

n  akk  J  (-  1),+J  det  (Hu  (i,  j))  j  >  i, 

k=i  / 


where  Hu  (r,  s )  is  as  before. 

Proof.  Denote  WU  by  E  =  (e,j).  It  is  clear  that  for  the  case  i  =  j,  E  =  In  where  In  is  the  nth  unit  matrix.  Now  consider  the 
case  j  >  i.  From  the  definitions  of  matrices  W  and  U , 

n  j 

eij  —  ^  ^  =  ^  ]  U)it&tj 

t=  1  t=i 


t=i+ 1  \k=i 
\2i+l  , 


= — + yi  ( n akk)  ab det (*> f)) 


_  ay  (-1)  l+  ai+hj  det  (Hu  (i,  i  +  1))  (- 1 )  I+  ai+2j  det  (Hu  ( i ,  i  +  2)) 

an 


aiiai+  l,i+l 


aiiai+1,i+1ai+2,i+2 


_  _  (— ^)i+J~2  aj-2j  det  (Hu  (i,j  —  2))  (-l)^"1  ahhj  det  (Hu  (i,j  -  f))  (_-l)i+J  det  (H„  (i,  j)) 


J— 2 

n  akk 

k—i 


J-l 

n  akk 

k=i 


J-l 

n  akk 

k=i 


=  n^)  (m-D2'  n  akk  +  ai+ ij(-l)2,+1  Yi  a'!k  det  (Hu  (i,  i  +  1)) 

\k=i  /  \  k=i+l  k=i+ 2 

j-l 

H - f  aj-2j  (-1),+J“2  Yi  det  0’>J  -  2))  +  (-1)I+J_1  aj-ij  det  (Hu  (i,  j  -  1)) 

k=j-\ 


+  (—  l)i+J  det  (Hu  ( ij )) 
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((  I!  a'*j  ’  Qi+t-ij  det  (Hu  (; i ,  i+t-  l))j  +  (-1),+J  det  (Hu  (i,j)) 

which,  by  Lemma  3,  implies  e,j  =  0,  and  the  proof  is  completed.  □ 

All  the  results  of  this  section  hold  for  lower  triangular  matrices,  with  the  obvious  modification  using  the  lower  Hessenberg 
matrix. 

Now  we  mention  an  interesting  fact  that  the  numbers  of  summed  or  subtracted  terms  in  computing  the  inverse  of  a  term 
of  an  upper  (lower)  triangular  matrix  are  the  generalized  order-/;  Fibonacci  numbers  defined  by 

fn=J2Cifn-i’  ”>0’ 

1=1 

where  /,_k  =  1  ,f*k  =  ■■■=/*  =  0. 

When  /;  =  2  and  Ci  =  c2  =  1,  the  generalized  order-2  Fibonacci  numbers  are  the  usual  Fibonacci  numbers,  that  is, 
f*  =  Fm  (mth  Fibonacci  number).  When  also  k  =  3,  Ci  =  c2  =  c3  =  1,  then  the  generalized  order-3  Fibonacci  numbers  by 
the  initial  conditions/^  =  1  ,/03  =  //,  =  0,  are 

1,  1,2,4,  7,  13,  24,  ..., 


which  are  also  known  as  tribonacci  numbers. 

Let  Un  =  (ify)  be  an  upper  (with  k  super-diagonals)  triangular  matrix  of  order  n  with  Ua  =  h,  for  1  <  i  <  n,  u,  ,+r  = 
Ei Li  Qr,i  for  1  <  i  <  n  —  r,  1  <  r  <  k  and  /i,’s  are  all  distinct  from  zero,  and,  arl' s  are  arbitrary. 

For  computing  the  inverse  of  L/„,  we  need  the  corresponding  Hessenberg  submatrices  defined  as  before  by  H  (r,  s). 
Therefore  we  should  note  that  the  numbers  of  summed  or  subtracted  terms  in  computing  the  inverse  of  a  term  of  Un  omitting 
the  signs  and  denominators  of  the  terms,  that  is,  the  number  of  required  summations  in  the  expansion  of  det  (H  (r,  s)), 
are  the  generalized  order-(s  —  r)  Fibonacci  numbers,  /ns_r.  To  show  that,  we  consider  a  (k  x  k)  upper  Hessenberg  matrix 

Hit  =  (a,j)  with  r-superdiagonals  whose  entries  are  given  by  a,-  +u  =  e,  for  1  <  i  <  k—  1,  au+r  =  Etli  K'!  for  1  <  i  <  k— r 
and  0  <  r  <  k,  et  ^  0  for  all  i.  If  one  superdiagonal  has  a  0  entry,  then  all  the  entries  in  this  superdiagonal  are  zeros.  That  is, 
if  al  i+r  =  0  for  some  i  and  r,  then  cr  =  0. 

Here  En  denotes  the  number  of  summed  or  subtracted  terms  in  detF/n.  For  example,  if 


H3 


a  b  c~ 
d  e  f 
.0  g  h_ 


then  detf/3  =  ahe  —  afg  —  bdh  +  cdg  and  so,  the  corresponding  E3  =  4. 

By  expanding  det  Hn  with  respect  to  the  first  row  without  any  simplification  in  the  entries  h" 1  =  EtLi  Qr,i.  we  get 

En  =  CiEji-1  +  c2En-2  +  •  •  •  +  cr£n_r. 


One  easily  computes  that  Ei  =  1,  £2  =  2,  E3  =  4,  . . . ,  Er  =  2r  ’. 

Consequently  one  can  see  that  £„  =  /nk+1  where /nk  is  the  generalized  order-/;  Fibonacci  numbers. 

For  example,  for  /;  =  4,  r  =  2,  let  Ci  =  2  and  c2  =  1,  that  is,  a,-  ,+i  =  Et=i  K'i  =  h'/i  +  k2 af,,+2  =  Et=i  = 
and  so 


H4 


>,(i) 


'  h'  ’ 
1.1  '  "2,1 

ei 

0 

0 


h(2) 

"l,1 

+h(V 

"l,2  '  "2,2 

e2 

0 


h (2> 

"1,2 

hw+h(V 

"l,3  '  n2,3 

e3 


0 

0 

hfl 

hlA  +  h 


(1) 

2,4_ 


Since  r  =  2,  Ci  =  2,  c2  =  1,  the  counting  sequence  £„  satisfies 


En  —  2En_i  +  £n_2, 


with  £  1  =  1,  E2  =  2,  so  En  =  Pn,  the  well  known  Pell  sequence. 

Therefore,  the  number  of  summed  or  subtracted  terms  while  computing  detH4  is  the  5th  Pell  number: 


det 


'/,(’)  +  /,(D 

"l.l  T-  l‘2A 

ei 

0 

0 

„  n  1,2  u2 


hn) 

"1.2 


ht2) 

,fl  1 

0 

0 

+  /i(1) 

^  "2,2 

h(2) 

"1,2 

0 

02 

h(1)  +h(1) 
"i,3  +  "2,3 

h(2) 

"i,3 

0 

e3 

h(i!l  +  f 

0) 
2,4 _ | 

1  1 , 2  /i  1 . 3 

kl,4  +  fll,l/l 

1.2^1,3^12,4 

+  /1 
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+  hl.2/l2. 1^1, 3^1, 4  +  hi,  1/11,2^2,3^2,4  ^1,1^1,3^2,2^12,4  +  lll.lh2,2hi  4h2.3  +  hi,2h2 ,lhlj3h2,4 

+  hi.2h21hi,4h2,3  +  h2,ihi,3h2,2hi,4  +  hi.1h2,2h2,3h2  4  +  hi,2h2.i h23h24  +  h21hi  3h2  2h2,4 
+  h2,ih2  2hi,4h2,3  +  h2,ih2,2h2  3h2,4  —  e4h2 4  —  e2h3jh2  2h14  —  e3hi,ihi,2h^  3 

—  eR  1hi,3h2,4  —  eR  lhi,4h2,3  —  e2h2  2h2jhj  4  —  e2h\  \h\  2h2  4  —  e3h4  4h2  3h2  2 

—  e3h\t2h2  4h\  3  —  e\h\  1h2,3h2,4  —  e2h2  2h2jh2  4  —  e3h2  jh2  3h2  2. 

From  the  above  example,  it  is  seen  that  there  are  29  terms  in  the  expansion  of  detLLi  which  is  the  5th  Pell  number,  P5. 


4.  The  inverse  of  an  r-banded  matrix 


In  this  section,  we  give  a  closed  formula  for  the  inverse  of  an  r-banded  matrix.  First,  we  consider  the  inverse  of  the  upper 
triangular  matrix  L\.  Flere  we  recall  a  well  known  fact  that  the  inverse  of  an  upper  triangular  matrix  is  also  upper  triangular. 
We  shall  give  the  following  lemmas  whose  proofs  are  straightforward. 


Lemma  5.  Let  the  lower  triangular  matrix  L  be  as  in  (4).  Let  £  =  (e,j)  denote  the  inverse  of  L.  Then 

e-l1’  .. 

,J  |  (— 1)I+J  det  (Le  (i,j)) ,  if  j  <  i, 
where  U  (i,  j)  is  defined  as  before. 


Lemma  6.  Let  the  upper  triangular  matrix  U  be  as  in  (5).  Let  G  =  (g,j)  denote  the  inverse  of  U.  Then 


Sij  = 


1 

w 

(— 1)I+J  det  (Lfu  (i,  j)) 

rR1 


if  j  =  i, 
if  j  >  U 


where  Uu  ( i,j )  is  defined  as  before. 

The  main  result  of  this  section  follows  from  the  LU  factorization  of  Br  n  and  the  previous  lemmas. 


Theorem  7.  Let  Dn  =  (dy)  denote  the  inverse  of  the  matrix  Brn 

lemmas.  Precisely,  if  we  let  S(i,  t,j)  =  det<u"(|JL)dePG(tj))>  f^en 

nr=i  kr 


.  Then  d, j  =  YTl=  \  where  git,  etJ-  are  defined  in  the  previous 


d 


ij  ~ 


,  1)l-+,-det(I/u(f,j)) 

n  a1 


1 

U} 

,vi 


n 

+  ^S(i,  t,  i), 

t>i 


+  <  — l)i+JT>(i-tj), 

t>] 


l(-l)i+Jdet (lt(i,j))  +  (-l)i+J^A(i,  t.j), 

t>i 


if  i  <  j 

if  i  =  j 
if  i  >  f 


Proof.  Since  B,„  =  LL1,  then  by  the  previous  two  lemmas,  we  get  D  =  U  !L  \  and  so, 

n  n 

d]j  =  /  ]  Sitetj  =  ^  '  Sitetj- 

t=l  t=max{i,j} 

By  taking  the  three  cases  i  <  j,  i  =  j,  i  >  j  and  using  the  expressions  of  glt,  etj  from  Lemmas  5  and  6,  the  claim  follows.  □ 


5.  Some  particular  cases 

If  we  take  r  =  2,  and  S2  n  to  be  a  symmetric  matrix,  and  label  a?  =  at,  b2  =  b~ 2  =  —  bi+1,  we  obtain  the  formulas  of  [11] 
for  the  inverse  of 
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■  Qi 

-b2 

o  - 

-fa2 

a2 

-fa3  ■  ■  • 

0 

0 

-fa3 

a3  -b4 

0 

S2,n  — 

0 

-bn- 

1  On—  1 

~b„ 

_  0 

0 

bn 

a„  - 

namely,  the  existence 

of  two 

sequences 

Uj,  Vj  such  that 

UtV  i 

0\V2 

UjV3 

UlVn 

UW2 

U2V2 

u2v3  ■  ■  ■ 

U2Vn 

s-1  — 

2,n  — 

u  tV3 

U2V3 

u3v3  ■  ■  ■ 

U3Vn 

_u,vn 

U2V„ 

u3vn  ■  ■  ■ 

The  sequences  u,,  Vj  can  be  determined  (as  Meurant  did  in  [11])  using  the  LU  decomposition  of  Theorem  1,  and  they  will 
depend  on  our  kj.,  mj.  We  will  not  repeat  the  argument  here. 

Consider  the  general  binary  recurrence  Gn+ 1  =  aGn  +  b  Gn_ j,  G0  =  0,  Gi  =  1.  Let 


-  1  0  . 

-a  1  0  ■  ■ 

—fa  —a  1  0 


°1 

0 

0 


0 

_  0 

Its  inverse  is  [12] 

T  Gl 
C2 
C3 


-a  1  0 

—fa  —a  1_ 


0  .  0  -| 

Ci  0  .  0 

G2  Gi  0  0 


Gn-i  Gn_ 2 
-  Gn  Gn_  i 


C2 

g3 


Gi  0 

C2  Gi_ 


In  the  same  manner,  by  going  through  our  argument,  one  can  obtain  many  other  known  (and  possibly  unknown)  matrix 
inverses  and  determinants.  We  end  by  displaying,  yet  another  example,  based  on  the  Chebyshev  polynomials  of  the 
second  kind,  Un(x)  =  ;  cos g  =  x,  which  satisfy  the  recurrence:  l/n+i(x)  =  2 xUn(x)  —  C„_i(x),  U0(x)  = 


1,  Lii(x)  =  2x.  The  symmetric  Toeplitz  S  = 


(-D 


j+jl  ttf— l  (a/2b)Un~jla/2b) 
U„(a/2b) 


has  the  inverse  [13]  S  1 i  =  (ty),j  where  %  = 


i  <) 


.  „(+,l  P/'-l  (a/2b)Un~i(a/2b)  .  . 

b  Un(a/2b ) 

It  could  be  interesting  to  use  a  variation  of  the  methods  of  this  paper  to  investigate  the  spectrum  of  general  r-banded 
matrices  (see  [14],  for  example). 
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